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Abstract 



D ' Massless quantum electrodynamics is studied at high temperature and zero chem- 

ical potential. We compute the Debye screening mass to order e 4 and the free energy 
s * , to order e 5 by an effective field theory approach, recently developed by Braaten 

and Nieto. Our results are in agreement with calculations done in resummed 
CN| ■ perturbation theory. This method makes it possible to separate contributions 

to the free energy from different momentum scales (order T and eT) and 
\ provides an economical alternative to computations in the full theory which 

involves the dressing of internal propagators. 
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1 Introduction 



Effective field theories both at zero and finite temperature have become of increasing 
interest in recent years [1-3]. Modern developments in renormalization theory have 
given meaning to non-renormalizable quantum field theories ("effective theories") 
and it is no longer obvious that renormalizabilty is an essential property of a use- 
ful field theory Q. The distinction between effective theories and "fundamental" 
(renormalizable) field theories becomes rather vague. Indeed, the modern view is to 
consider every field theory as effective, valid and with predictive power at a certain 
scale. 

There has been a tremendous progress in perturbative calculations of the ther- 
mal properties for quantum field theories at high temperatures since the pioneering 
work of Dolan and Jackiw some twenty years ago |J. Today, one can essentially 
distinguish between two methods for computing high temperature properties of a 
quantum field theory. The first method is based on resummed perturbation the- 
ory. Resummation is a reorganization of the perturbation series in order to be 
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actually perturbative. Resummation gives rise to effective thermal masses (and in 
non-Abelian theories effective vertices, too) that provide an infrared cut off in the 
loops. The dressed propagators are thus well behaved and the use of these in loops 
yields the contributions to the free energy that are non-analytic in g 2 , where g is the 
coupling constant. The use of resummed perturbation theory amounts to summing 
an infinite set of higher order diagrams. This resummation scheme has largely been 
developed by Braaten and Pisarski |6| and has been applied to QCD by a number 
of authors in connection with the gauge dependence of the gluon damping rate and 
the screening mass (Refs. [7,8] and references therein). The method has also been 



used to to calculate the free energy in QED [9,10] and in the standard model fll 



The second method is based on dimensional reduction [|l]|. The crucial obser- 
vation is that the only contribution to a correlator that is not exponentially damped 
at scales larger than T comes from the static bosonic modes. One strategy is there- 
fore to integrate out the non-static bosonic modes as well as the fermionic modes 
to construct an effective three-dimensional field theory of the zero modes, in which 
the masses and couplings are temperature dependent. This approach has been used 
in the study of high temperature QCD and QED [12| as well as in investigations 
of phase transitions in spontaneously broken gauge theories [1,13] (in particular the 
electro- weak phase transition). 

Instead of explicitly integrating over the heavy modes (which produces non- 
local terms beyond one- loop [0), one simply writes down the most general three- 
dimensional Lagrangian, that respects the symmetries of the system. One then 
computes static correlators in the two theories and require that they match. This 
matching requirement can actually be taken as the definition of the effective field 
theory and provides the relationship between the coupling constants in the effective 
theory and the underlying theory. This method has previously been applied to 
massless 4 theory |2| and QCD [15|]. In the latter case it does not only provide 



a convenient way of calculating the free energy in the high temperature limit, it 
does also solve the infrared catastrophe of Non-Abelian gauge theories. It is a well 
known fact that the free energy of Non-Abelian gauge theories may be calculated 
to fifth order in the coupling using resummed perturbation theory []16 |. However, 



the method breaks down at order g 6 and this has been interpreted as a sign that 
the transverse gluons have a mass of order g 2 T, which cannot be calculated in 



perturbation theory JTTJ| . This magnetic screening mass should then provide the 



necessary infrared cutoff. 

In the present paper we would like to apply the method to QED and confirm 
some recent results obtained by calculations in the full theory [9,10]. Our work will 
serve as a nice demonstration of the efficiency of the effective field theory approach 
by allowing one to work with a single scale at a time. 
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The outline of this work is as follows: In section two we compute the screening 
mass to fourth order and the free energy to fifth order in the coupling constant 
in high temperature QED. In section three we summarize. Throughout the work 
we use the imaginary time formalism and the theory is regularized by working in 
4 — 2e dimensions together with the MS regularization scheme. We use the following 
shorthand notation for the sum-integrals that appear below: 

J{dK} f (K) = J4^m), (1) 



fc =27r(n+l/2)T 



where K 2 = k\ + k 2 . In the Feynman diagrams a solid lines denotes a fermion, a 
wavy line a photon and a dashed line a ghost. 



2 QED at High Temperature 



The partition function can be written as a path integral 



dr / d xC 



Z = J VAp Vip Vip exp 
where the Lagrangian of QED reads 

£ = -F^F^ + ^7"(0„ - ieA^jip + £ gf . 



(2) 



(3) 



In this paper we choose to work in Feynman gauge where L g f = ^(d^A^) 2 . In the 
effective theory the partition function reads 



Z = e - f{A)v [ VAiVA exp \ - f d 3 x£ eS 



(4) 



Here C e Q is the effective three-dimensional Lagrangian to which we shall return 
shortly. The prefactor /(A) can be interpreted as the coefficient of the unit operator 
in the effective field theory. It depends on the ultraviolet cutoff A (which is intro- 
duced to regularize the effective three-dimensional theory) in order to cancel the 
A dependence in the path integral in Eq. (^j) 0. The effective three-dimensional 
theory consists of a gauge field and a real massive self interacting scalar field in the 
adjoint representation of the gauge group: 



C 



1 • • ' •■ \2 , ' ...2, \ , 2 



F l3 F l3 + ^{d lP f + ^m 2 (A)p z + A(A)p 4 + C gi + 5C. 



(5) 



At this point a few comments are in order. Firstly, we remark that the scalar field p 
is identified with the time component Aq of the gauge field in the full theory, up to 
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normalizations. This implies that the field Aq does not interact with the gauge field 
in the effective theory, since QED is an Abelian gauge theory. The fact that A may 
develop a thermal mass is simply a consequence of the lack of Lorentz invariance at 
non-zero temperature. Secondly, it can easily be demonstrated from the matching 
of four-point functions that A(A) is of order e 4 and that the term A(A)p 4 does not 
contribute to the Debye mass or the free energy to the order we are calculating (This 
interaction contributes at order e 5 to the screening mass and at order e 6 to the free 
energy) . 

Furthermore, 5C represents all local terms which respect the symmetries of the 
theory such as gauge invariance. This includes renormalizable terms, such as g(A)p e , 
as well as non-renormalizable ones. 



The mass parameter. In order to obtain the free energy to order e 5 , we need to 
compute the short distance coefficient m 2 (A) to order e 4 . The e 3 and e 5 contri- 
butions to In Z is namely given by a one-loop calculation in the effective theory 
(see below). The simplest way to determine m 2 (A) in the effective theory is by 
demanding that the screening mass in the two theories match. The screening mass 
of the particles is defined as the location of the pole of the propagator for spacelike 
momentum [0: 

p 2 + n 00 (0,p) = 0, p 2 = -ml (6) 

This definition is gauge fixing independent, which one can prove on an algebraic 
level |17| . This is in contrast with the one found in the book by Kapusta |TB[ and 
which commonly used in the literature!: 

m 2 = n 00 (p = 0,p^0). (7) 

Normally, as in the present case, the two definitions are equivalent to leading order in 
the coupling constant. Beyond leading order they do not coincide. The requirement 
above implies that 



p 2 + m 2 (A) + n eff (p, A) = 0, p 2 = -m 2 , 



where U e g (p, A) is the self-energy of p in the effective theory. We shall do the 
matching in strict perturbation theory. Generally, this means that the expression for 
the screening mass differs from the one obtained in resummed perturbation theory, 
which correctly incorporates the effects of electrostatic screening. The solution m 2 
is of order e 2 , which implies that one can expand II (p 2 ) = IToo(0,p) in a Taylor series 
around p 2 = 0. To determine the screening mass (in strict perturbation theory) to 



1 In quantum electrodynamics the latter is gauge fixing independent, but it is not renormalization 
group invariant. In non- Abelian theories it is also gauge dependent. This is to be expected since 
it is an off-shell quantity Q. 
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order e 4 , we must calculate IT(p 2 ) to one loop order and IT(0) to two loop order, 
and the screening mass is then given by 



11(0) l-n'(O) 



(9) 



The symbol ~ indicates that Eq. (|9|) only holds in strict perturbation theory. Now, 
strict perturbation theory in the effective theory means that the mass term should 
be treated as a perturbation. The corresponding contribution to II e g- is shown in 
Fig. [I], where the blob indicates the mass insertion. The matching relation simply 
becomes m 2 (A) « m 2 . 



The self-energy to one-loop order in the full theory reads 



n(p 2 ) = e 2 J{dK}Tr{ l0 ^ l0 



{dK}^ + e* 4 -pi 



(10) 



The corresponding Feynman diagram is shown in Fig. |2|. The sum-integrals in Eq. 
( |10|) are standard and can be found in e.g. Ref. [jn|. Carrying out the renormaliza- 
tion of the wave function in the usual way, one finds: 



n( P 2 



e 2 T 2 



+ -^(27^-1 + 2^-^-; 

127T AlTl 



'111 



Here A is the scale introduced by dimensional regularization. The two-loop expres- 
sion for the self energy at zero external momentum can be found either by a direct 
computation of the two- loop graphs or by applying the formula (see Ref. |T8| ) 



n(o) 



, o 2 p 

<9/i 2 ' 



(12) 



where P is the pressure and \i is the chemical potential. This requires the calculation 
of the free energy to two-loop order including the chemical potential. This can be 
carried out by contour integration [18|. Using the result —T 2 e A /8n 2 for the e 4 
contribution, we find 



e 2 T 2 e 2 v 2 
n(p 2 ) = ^ + ^(2 7 *-l + 21n 

The mass parameter to order e 4 then becomes: 



A , 
4ttT' 



m 2 (A) = T 2 



e 



e 4 A 
(2 TB -l + 21n 



367T 



AttT' 



8tt 2 



(13) 



2^4 



T 2 e 
~8~7 2 



(14) 
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By using the renormalization group equation for the coupling constant, 

one can easily demonstrate that Eq. fll4]) is independent of A. Thus, up to correc- 
tions of order e 6 , we can replace A by an arbitrary renormalization scale fi. 

The coupling constant. We would also like to make a few remarks about the gauge 
coupling e% in the effective theory. Using the relation between the gauge fields in 
the two theories 

4 3 = VfA l} (16) 

one finds that e| = Te 2 to leading order in e 2 . There is no dependence of A at this 
order. Beyond leading order one must compute and match four point correlation 
functions 0. 

The coefficient of the unit operator. We shall now compute /(A) to order e 4 in 
strict perturbation theory. A strict perturbative expansion in the coupling constant 
is normally plagued by severe infrared divergences, which are caused by long-range 
forces. Physically, these forces are screened at the length scale l/(eT) and can be 
taken into account only by reorganizing perturbation theory which amounts to sum- 
ming an infinite set of loops. Nevertheless, the strict perturbation theory can be 
used to compute the short distance coefficient /(A) and we shall do so by matching 
calculations of In Z in the full theory and in the effective theory. From eqs. (|2|) and 
(f|), we see that the matching condition reads 

lnZ = -/(A)F + lnZ eff . (17) 

The calculation of InZ in the full theory involves one-loop, two-loop and three-loop 
diagrams. The graphs are displayed in Figs. [3]— The diagrams can be computed 
by the methods developed in Refs. [9,19] and details may be found there. We find: 

T\nZ 7T 2 T 4 r ll 5e 2 e 4 / 20 , , A , 8C'(-3) 16 C'(— 1) 

+ 7^— d --ln(-— ) + - S 



V 9 



20 32tt 2 256tt 4 V 3 HnT J 3 C(-3) 3 ((-1) 

at 



319 208, n 

- 4>y E 1 In 2 

' 12 5 



Here ~ again indicates that Eq. ( |T8"D only holds in strict perturbation theory, and 
A is the momentum scale which is introduced by dimensional regularization. A may 
again be traded for an arbitrary \i by RG-arguments. 

We now turn to the effective theory. The mass parameter is viewed as a pertur- 
bation in the effective theory. lnZ e g is then given by one-loop contributions from 
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the gauge field and scalar field (Figs. |3|— ||), plus an additional one-loop diagram 
with a mass insertion (which is indicated by a blob in Fig. ||). The computation 
is rather simple since one-loop contributions involving massless fields vanish iden- 
tically in dimensional regularization (see Eq. (|2lD below) and so does \nZ e Q. The 
matching condition therefore turns out to be 

^ - -f(m (19) 

and /(A) is given by the right hand side of Eq. (|TE|). The function F = f(A)T can 
be viewed as the contribution to the free energy from the short distance scale 1/T. 
With the comments after Eq. (fT8| ) in mind, it is clear that /(A) has no dependence 
of A at the order we are calculating. 

Now that we have determined the short distance coefficients m 2 (A) and /(A) in 
the effective theory, we can calculate the screening mass as well as the free energy 
in QED to order e 4 and e 5 , respectively. 

The screening mass. In order to calculate the physical screening mass, one must 
include the mass parameter m 2 (A) in the free part of the Lagrangian. Furthermore, 
since the term A (A) goes like e 4 , the 4 term does not contribute to the screening 
mass to order e 4 , as we mentioned above. This implies that the physical screening 
mass! is equal to the short distance coefficient m 2 (A) to order e 4 : 

re 2 p 4 u T 2 P 4 i 

m 2 s = T 2 [\- ^- 2 (2 lE - 1 + 2 In JL) - Lt. . (20) 

It is easily checked that the result is RG-invariant as required. Furthermore, our 
result agrees with the calculation of Blaizot et al to order e 4 |pO| . Note also that 



there is no e 3 term in the expression for the screening mass in contrast with both 
4 theory and SQED. The reason is that there are no bosonic propagators in the 
one-loop self-energy graph in QED and fermions need no resummation, since their 
Matsubara frequencies are never zero. 

The free energy. The calculation of the free energy in the effective theory is now 
straightforward, it is simply a one-loop computation. However, we must now include 
the physical effect of screening, which amounts to consider the mass parameter as 
a term in the free part of the Lagrangian. Using dimensional regularization (see 
e.g. Ref. ET[) the one-loop integrals are perfectly finite after regularization and 



2 It is also equal to the screening mass obtained in the strict perturbation theory to order e 4 . 
This is due to the fact that there is no e 3 term, as explained above. At order e 5 and higher they 
differ, caused by the dressing of bosonic propagators in the two-loop diagrams. 
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independent of the renormalization scale: 



rl 3 h 1 

\n(k 2 + M 2 ) = — — -M 3 . (21) 



(2tt) 3 v ' 6tt' 

The contribution from the gauge field vanishes. Using the expression for the mass 
of the scalar field and expanding it in powers of e yields the following contribution 
to the free energy: 



\nZ eS = T 4 



3 5 

— e -= ^^(41n-^ + 4 7i; + 81n2 + 9 

36^ 7T 576^3tt 3 V 4ttT 



:22i 



This term takes into account the effects from long distance scales of order l/(eT). 
Note that it is non-analytic in e 2 . Furthermore, in resummed perturbation theory 
the e 3 and the e 5 terms would arise from the dressing of the photon propagator in 
the two and three-loop diagrams, respectively. Putting the results from Eq.s (|18D 
and (|22|) together, one obtains 

IItt 2 5e 2 e 3 e 4 / 20 /i 8 C'(-3) 16 C'(-l) 



InZ = T 



180 248 3QV3ti 2304tt 2 V 3 v 4ttT 7 3 C(-3) 3 CH 
_ 47E _^ + ^ ln2) __fL_ (4lD _£_ +47 , + Sl„2 + 9 )]. 

This is the main result of the present paper and it is in agreement with the compu- 
tation of Zhai and Kastening |J, and Parwani |IIJ] who use resummed perturbation 
theory. Eq. (23) is correct as well as renormalization group invariant up to cor- 
rections of order e 6 In e. The latter property may easily be checked by using the 
one-loop /3-function in QED. 

Finally, we would like to comment upon a computational difference between 
QED and QCD. In QCD the computation of the free energy involves the construction 
of two effective field theories, which reflects the fact that there are contributions from 
three different momentum scales (T, gT and g 2 T, where g is the gauge coupling) 



T5j. The first effective field theory, called electrostatic QCD (EQCD), consists of 
the magnetostatic A? field and the electrostatic field Aq. The unit operator /qcd as 
well as the mass parameter in EQCD is then determined by the matching procedure 
and /qcd gives the contribution to the free energy from the short distance scale 1/T. 
The second effective field theory is called magnetostatic QCD (MQCD) and consists 
simply of the self interacting magnetostatic gauge field Af. Again the unit operator 
of this effective theory, /eQCD> can be determined and yields the contribution to 
the free energy from the distance scale l/(gT). Now, the perturbative expansion in 
MQCD is plagued with infrared divergences, implying that the functional integral 
can only be calculated non-perturbatively, e.g. by putting MQCD on a lattice. 
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Using lattice simulations the path integral may be computed, so that one obtains 
the contribution to the free energy from the scale l/(g 2 T). One can, of course, also 
construct a second effective field theory in QED, but it is completely unnecessary. 
The restriction to the magnetostatic zero modes simply yields a free photon field 
theory in three dimensions. This in turn implies that there are no contribution from 
momentum scales e 2 T in QED. The same remark applies to scalar electrodynamics. 



3 Summary 



In this work we have computed the free energy at high temperature in massless 
quantum electrodynamics to order e 5 using an effective field theory approach. This 
approach makes it possible to separate contributions which come from different 
length scales (1/T and l/(eT)). Furthermore, this method involves considerably 
less effort in obtaining the free energy to a given order in perturbation theory than 
computations in the full four dimensional theory. 



The author would like to thank F. Ravndal for comments and suggestions, and 
H. Haugerud for carefully reading the manuscript. 
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FIGURE CAPTIONS: 

Figure 1: Self-energy correction in the effective theory. 
Figure 2: One- loop self-energy correction. 
Figure 3: One-loop diagrams. 
Figure 4: Two-loop diagram. 
Figure 5: Three- loop diagrams. 
Figure 6: One- loop diagram with mass insertion. 
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